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works, the extracted inelastic overlap function presents a peripheral effect (tail) above 2 fm and the 
extracted opacity function is characterized by a zero (change of sign) in the momentum transfer 
space, confirming results from previous model-independent analyses. Analytical parametrization 
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posed representation is commented. A short critical review on the use of Martin's formula is also 
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I. INTRODUCTION 



The theoretical description of the high-energy elastic hadron scattering still constitutes a hard 
challenge for QCD: perturbative techniques cannot be formally applied at this soft (large distances) 
region and nonperturbative approaches are not yet able to predict soft scattering states from first 
principles without model assumptions P, Q some of them unjustified. On the other hand, in the 
phenomenological context, a wide variety of models can describe the experimental data, but since 
they are characterized by different physical pictures Q a comprehensive, well accepted and unified 
formalism is still missing. This situation has to be contrasted with the great expectations from the 
new experiments at the CERN LHC, in special the investigation of the soft diffractive processes in 
proton-proton (pp) collisions by the TOTEM Collaboration 0, Q . 
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At this stage, model-independent analyses, aimed to extract from the experimental data empirical 
information on what is theoretically unknown, or on what cannot yet be evaluated, may be an impor- 
tant strategy for realistic and necessary developments of the theory in the soft sector. As first steps 
in that direction, we have investigated the elastic scattering based only on some general principles, 
high-energy theorems and model-independent fits to the differential cross section data [6l— Il3jj . Since 
the impact parameter and the unitarized eikonal representations still constitutes a useful framework 
for model developments, the aims have been to extract empirical information on the overlap and 
eikonal functions (impact parameter and momentum transfer spaces). However, since these analyses 
demand, as input, experimental information on the differential cross section at intermediate and large 
values of the momentum transfer, they have been limited to proton-proton scattering in the energy 
interval 19.4 - 62.5 GeV (as obtained from the CERN SPS, Fermilab and CERN ISR, nearly forty 
years ago). Therefore, the new data from the LHC on pp collisions at 7 - 14 TeV, covering the region 
up to large momentum transfer 0, Q will certainly be crucial for further developments. 

One of the main problems with these model-independent analyses is the fact that the phase of 
the scattering amplitude constitutes a physical observable only in the forward direction (at zero 
momentum transfer), namely we do not have information on the contributions from the real and 
imaginary parts of the amplitude beyond the forward direction. As a consequence, if we look for 
empirical analyses, unbiased by any model/phenomenological ideas @ different kinds of contributions 
must be tested and reliable information should not depend, or not strongly depend on the details of 
these contributions. With this strategy in mind, we have already discussed two different empirical 
solutions for data reductions both statistically consistent For example, the constrained results 

obtained in [9J are in agreement with the standart picture, characterized by the dominance of the 
imaginary part of the amplitude, except at the dip position, where the imaginary part presents a zero 
(change of sign) and the dip is filled up by the real part. On the other hand, the recent unconstrained 
results of @ indicate the dominance of the real part at intermediate values of the momentum transfer, 
just above the dip region. 

Motivated by these different results, both consistent on statistical grounds, we consider here a third 
possibility with different contributions from those obtained in the above mentioned previous analyses 
and developed by means of a novel almost model-independent representation for what has been know 
as Martin's real part formula This formula connects real and imaginary parts of the elastic 

scattering amplitude in terms of the transferred momentum in the collision; it has been derived in the 
context of high-energy theorems and is embodied by a scaling property. 

In this work a representation for this formula, without the full scaling property, is used in model- 
independent analysis of elastic pp scattering at high energies. By means of an empirical parametriza- 
tion for the imaginary part of the amplitude, the use of the Martin's formula for the evaluation of the 
corresponding real part, and fits to pp differential cross section data (in the interval 19.4 - 62.5 GeV), 
the overlaps (elastic, inelastic and total) and eikonal functions are extracted in an almost model- 
independent way. The reproduction of all the experimental data is quite good and the two main 
results are: the extracted inelastic overlap function (impact parameter space) is characterized by a 
peripheral effect (tail) above 2 fm and the extracted imaginary part of the eikonal (opacity function) 
in the momentum transfer space presents a zero (change of sign), in the region 6-7 GeV 2 . Analytical 
parametrization for these empirical results are introduced and discussed in certain detail. 

In the context of our global investigation on the inverse problems in elastic hadron scattering 
, Q the focus here is only in the foundations, applicability and some consequences of the proposed 
representation. With this aim we shall only outline references and discussions related to results from 
previous analyses and phenomenological models. A global critical review, taking into account results 
that are common to all the analyses (references @, Q and those presented here) , as well as detailed 
discussions on the implications in the phenomenological context will be presented in a forthcoming 
work. 

The manuscript is organized as follows. In Sect. II we review the original derivation of the Martin's 
formula, present some critical comments on subsequent developments and introduce our representa- 
tion. In Sect. Ill we describe the fit procedures in some detail and display the results of the data 
reductions. In Sect. IV we treat the extracted quantities, the inelastic overlap function and the 
eikonal in the momentum transfer space, making reference to results from previous analysis and some 
phenomenological models. The conclusions and final remarks are the contents of Sect. V. 
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II. MARTIN'S FORMULA 



General model-independent properties of the scattering amplitude can be obtained from fundamen- 
tals principles of local quantum field theory and a set of rigorous theorems [15l - ll8l |; in this context, 
unitarity, analyticity and crossing play a central role. Although rigorous, some results are obtained 
under specific conditions, which define and limit their range of validity and applicability. Martin's for- 
mula has been obtained in this context and therefore it is important to stress the conditions involved. 
In this Section we first treat the assumptions and main steps in the original derivation [bH . [l9j , followed 
by a short critical review on the subsequent interpretations and applicability of the formula. After 
that, we introduce a representation for almost model-independent analysis, with a detailed discussion 
on the conditions involved. 



A. Original derivation 



For elastic processes, to + to — > to + to, the laboratory energy E is a useful crossing symmetry 
variable. Let F be the scattering amplitude and t the square of the momentum transfer. In the 
high-energy limit (E >> to) the optical theorem reads [201 ] 

lmF(E, t = 0) 
hi 

and the Froissat-Martin bound states that [2l| - [23l | 

o-tot <C[\nE} 2 , 

with C a constant. 

In what follows we shall use the symbol ~ for an asymptotic equality to with a constant factor. 
Auberson, Kinoshita and Martin [24] (AKM) have shown that if the above bound is saturated (<T to t ~ 
[In E] 2 ), then for 

k 

\t\ < t; v7 , k arbitrary constant, (1) 

[In E\ 1 



an scaling property holds 



F(E,t) 



f(r), T = t[\nE} 2 (2) 



F(E,0) 

and therefore, in that case, from the optical theorem, 

F(E, t) - E[\nE} 2 f(t[\nE} 2 ). 

Real and imaginary parts of the amplitude can be correlated through high-energy theorems based 
on analyticity and the crossing relations established by Bros, Epstein and Glaser in the context of 
quantum field theory [25L [26j . The main ingredient is the condition of polynomial boundedness and 
the Phragmen-Lindeloff theorem which leads to the asymptotic uniqueness [27]]. By assuming that 
the odd-under-crossing amplitude can be neglected at high ener gies , exact crossing for the even (+) 
amplitude demands that the above real amplitude has structure [l9l l27j 



F + ~iE[lnE-i-] 2 f(t[lnE-i-} 2 ). 



7T 

Since [InE — i— ] 2 pa [lni?] 2 — i7rlni?, the function / can be expanded and in first order we get 



17TT ITTT df 

fiT ~^E )Kif{T) -^E^' 
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so that in the high-energy limit the Martin's result is obtained 

JmF + (E,t) ~ E{\nE} 2 f(T), (3) 
ReF + (E,t) ~ 7rE\nE^-(Tf(T)). (4) 

Hence, for a given /(r) (imaginary part), the real part can be evaluated. We stress the four main con- 
ditions in this derivation: (1) asymptotic limit (Froissart-Martin bound reached); (2) even amplitude; 
(3) scaling property; (4) \t\ < fc/[ln£;] 2 . 



B. Geometrical Scaling, Developments and Critical Comments 



"Sometimes these theorems called 'asymptotic' are not expected to apply to present avail- 
able energies, but in fact some of them do apply already now. Then they are no longer 
theorems but they give useful trends." 

A. Martin, G. Matthiae [l]| 

In what follows we shall consider the center-of-mass system, with focus in this section on pp and 
pp elastic scattering. Let A(s, q) be the complex amplitude, where s is the square of the energy and 
q 2 = —t. The physical quantities of interest here, with the corresponding normalizations, are the 
differential cross section 

^=<x\A(s,q)\ 2 , (5) 

the total cross section 

otot = 4irTmA(s, 0), (6) 



and the p parameter, 



ReA(s,0) 

p{s) = wm) ■ (7) 



Martin's result (3-4) was published in 1973, a few months after the introduction of the Geometrical 
Scaling (GS) by Dias de Deus [2^|. This phenomenological approach, further developed by Dias de 
Deus, Buras and Kroll [29l - [3l| was based on the empirical evidence of a scaling in elastic pp scattering 
at the ISR energy region (rj 20 — 60 GeV) in terms of the variables 



1 da 



= $(^(7tot). (8) 



CTfot dq 2 

The geometrical aspect concerns the fact that for (T to t(s) oc R 2 (s), with R(s) an effective interaction 
radius, the increase of the total cross section and all the energy dependence involved in da/dq 2 is only 
due to an expansion (geometrical) effect of the hadron. 

Now, since from first order derivative dispersion relations [32T - [3^ | and for an even amplitude, 

P(s) ~ ~ ; otot(s), (9) 

2cr t ot(s) dins 

the Martin results, Eqs. (3) and (4), along with normalizations (5), (6) and (7), can be put in the 
form 



+ [/(r)] 2 j, (10) 

which also characterizes the GS model (29l - [3l| . Eq. (8), for r = q 2 otot (not necessarily asymptotic 
energies). 

Therefore, the result expressed by Eq. (10) (hereafter referred to as Martin's formula, as in the 
literature), has two independent foundations: 



da 
dq^ 



16tt 
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(1) the AKM scaling, Eq. (2), a rigorous formal result, deduced under strictly conditions of asymp- 
totic energies (Froissart-Martin bound saturated), valid in a limited interval of momentum trans- 
fer, Eq. (1) and for even amplitude; 

(2) the GS, an empirical result verified at the ISR energy region (certainly not asymptotic), without 
restriction in momentum transfer or crossing properties of the amplitude. 

That seems a peculiar and/or accidental fact, putting in evidence some intrinsic difficulties asso- 
ciated with the asymptotic concept and the elastic processes in hadronic interactions, in general, as 
commented by Martin and Matthiae [l9j . Let us shortly review some further developments concerning 
the energy and momentum transfer variables associated with the Martin's formula. 

After the introduction of these scaling in the seventies, the new energy domain reached with sub- 
sequents experiments on pp scattering at the CERN SppS collider (546 GeV) and Fermilab Tevatron 
(1.8 TeV), have shed light on some aspects involved, as follows. The GS predicts a constant value for 
the ratio between the elastic and total cross section, a result in agreement with the experimental data 
at the ISR energy region. However, at the Collider and Tevatron energies it has been verified that 
this ratio increases [331, leading to the breakdown of the GS. On the other hand, in the phenomeno- 
logical context, Henzi and Valin have shown that from the ISR to Collider, the energy evolution of 
the dip-bump structure of the differential cross section demands not only a geometrical expansion but 
also blackening and edging effects [36l.l37j. With their BEL approach (Blacker, Edgier, Larger), based 
on dispersive diffraction theory, corrections have been introduced in the Martin's formula, leading to 
extensions at higher energies [38j. 

Another aspect concerns the interval in momentum transfer in which Eq. (10) is expected or 
supposed to be valid and that plays a central role in our analysis. In this respect and in the context 
of the AKM scaling, Kundrat and Lokaji'cek have developed a detailed numerical analysis [39M4l1 | 
showing that the function /(t) in (10) is real only in a limited interval of the scaling variable r, 
corresponding to the region in momentum transfer < q 2 < 0.15 GeV 2 . Therefore, the use of the 
Martin's formula outside this region has no physical meaning. This conclusion, however, has been 
criticized by Kawasaki, Maehara and Yonezawa [42[ and according to these authors, Martin's formula 
can be directly deduced through derivative dispersion relations in the case that GS holds, which is 
verified, at least, at the ISR energy region. Therefore, in this energy domain there would be no limit 
for the applicability of the Martin's formula in terms of the momentum transfer. 

Despite the above discussion, a crucial result for our purposes has been introduced by S.M. Roy 
[IH in the investigation of unitarity inequalities connecting Re A(s,q) and Im A(s,q). The main 
ingredient, suggested by the author and raised by A. Martin (see acknowledgments in |43|), concerns 
the possibility of a non-scaling only in the imaginary part of the amplitude. In fact, if Re A is small 
compared with Im A it seems reasonable to consider the scaling property only in the evaluation of 
the real part. The proposal by Roy has been to consider Eq. (10) but now with the substitution 

ei \ i \ ImA(s,q) 

that is, the scaling is supposed to hold only in the small contribution of the real part. 

Although under limited formal condition of validity, this proposal can be assumed as a representation 
for the Martin's formula applicable, in principle, to all values of s and q 2 . In fact, despite its "hybrid" 
character the efficacy of this representation in the phenomenological context has been demonstrated, 
subsequently, by several authors, for both pp and pp scattering and different regions of the momentum 
transfer and energy with experimental data available [44T - I53T ] . 



C. A representation for almost model-independent analyses 

As commented in our introduction, motivated by two different model-independent results for the 
contributions from the real and imaginary parts of the amplitude @, Q we consider here a third 
possibility, based on a representation for the Martin's formula. The point is the Roy's proposal, Eq. 
(11), with a model- independent parametrization for the imaginary part of the amplitude, in the usual 
form of a sum of Gaussian in q (or exponential in q 2 ). Specifically, in Eq. (11) we consider 



n 

Eh 2 



(12) 
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where a.; and hi, i — 1, 2, n are real free parameters. With this assumption and through Eq. (10), 
we obtain the following representation for the full complex amplitude 



A(s,q) 



+ i 




(13) 



Using as input the experimental values of 0t o t and p in each energy analyzed, we shall consider fits 
to the differential cross section data through Eqs. (5) and (13). We note that in this representation s 
and q 2 are independent variables and the number of free parameters is 2n — 1, due to the constraint 



(14) 



However, before considering the above parametrization in our data analysis, some critical comments 
on its applicability here are in order: 

• The analysis by Kundrat and Lokajfcek [39 41] concerns the scaling variable r = q 2 In 2 s or 
r = q 2 (Jtot- Since in the representation (13) the variables s and q 2 are independent (Roy's 
proposal), we understand that the restriction to the narrow interval in the momentum transfer, 
referred to in Sect. II. B, does not apply. 

• As told before, we shall treat only pp scattering in the energy interval 19.4 GeV < ^/s < 62.5 
GeV, a region (ISR) where the GS is empirically verified [54j. Therefore, we see no restriction 
on its use in the evaluation of the real part of the amplitude. 

• On the other hand, just the above energy interval brings about some shortcomings related to 
formal and practical aspects involved. In fact, at the ISR region, pp and pp scattering are 
distinct in what concerns the experimental data on da/dq 2 , cr t ot and p. Therefore the neglecting 
of the odd-crossing amplitude (connected with the asymptotic limit in the Martin's derivation) 
cannot be justified. Now, even if the odd (-) contribution could be taken into account, that 
would demand a simultaneous analysis of pp and pp scattering, since 



A- 4- A 

_ ^PP ^ ^PP 



(15) 



However, in that case model-independent fits cannot be statistically developed due to the small 
interval in momentum transfer with data available on pp scattering [1 21 ] . 

Therefore, despite the advantages of parametrization (13), we understand that it has a limited formal 
justification in what concerns the evaluation of the real part of the amplitude for pp scattering in the 
energy interval of interest here. That was the reason why we refer to an almost model-independent 
analysis. 

On the other hand, we stress that we do not intend to extract any direct physical implication 
or interpretation from the specific contribution of the real part of the amplitude. On the contrary, 
as commented in the introductory section, our strategy is to look for global properties from different 
parametrization and fit results that are not connected with or do not strongly depend on the (unknown) 
contributions from Re^4 and Imj4 beyond the forward direction. With that in mind and as we shall 
show in the following sections, parametrization (13) brings new useful insights in the investigation of 
the inverse problems in elastic hadron scattering. 



III. FIT AND RESULTS 



Our purpose is to use parametrization (13) to fit differential cross section data through Eq. (5), 
using as input the experimental values on crtot and p in each set analyzed. In this section, we first 
refer to the data set in which our analysis is based (Sect. III. A), then we discuss the fit procedures in 
certain detail (Sect. III.B) and after that we present the fit results and some critical comments (Sect. 
III.G). 
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A. Experimental data 

In the context of our empirical analyses , the evaluation of the uncertainties in the extracted 

quantities plays a central role. Since the fits are based on the differential cross sections, that demands 
sets with available data covering the largest intervals in momentum transfer, mainly in what concerns 
the region of high momentum transfer. As already discussed in some detail [f| d, [l2j this condition 
limits the analysis to six sets of pp elastic scattering data, at yfs = 19.4 GeV (from the Fermilab and 
CERN SPS) and yfs = 23.5, 30.7, 44.7, 52.8 and 62.5 GeV (from the CERN ISR). In addition, we 
make use of the empirical result that at the ISR energy region (23.5 - 62.5 GeV) the differential cross 
section data above q 2 ~ 4 GeV 2 do not depend on the energy 55-58]. This fact allows the inclusion in 
each set (ISR) the data at 27.4 GeV (from Fermilab), covering the large momentum transfer region: 
5.5 < q 2 < 14.2 GeV 2 (see @ for a quantitative discussion on this respect). 

The data set have been collected from final published results and a complete list of references, with 
comments, can be found in Q for the data at 19.4 GeV and in Q for the ISR data (23.5 - 62.5 GeV). 
The differential cross sections data include the optical point, 



da 
dq^ 



-Ld + P 2 ) (16) 



and the data above the Coulomb-nuclear interference region, q 2 > 0.01 GeV 2 . All the available 
data have been included, without any kind of selection or exclusion of data points. As a test for 
goodness of fit we shall consider the \ 2 P er degrees of freedom (DOF) and since this test is based on 
the assumption of a Gaussian error distribution, only the statistical errors of the experimental data 
are taken into account (we do not include systematic errors). We shall return to this point in our 
comments at the end of Sect. III.C. 



B. Fit procedures 

The data reductions have been performed using the CERN-Minuit code [5!| through successive runs 
of the MIGRAD minimizer and with the confidence level for the uncertainties in the free parameters 
fixed at 70 % . The nonlinearity of the fit demands the start values of the free parameters and looking 
for a completely unbiased procedure the following methodology has been used. 

We have begun the analysis with the differential cross section set corresponding to the largest 
experimental information available, namely y/s — 52.8 GeV (plus data at 27.4 GeV, as explained in 
Sect. III. A) and through the following steps: 

1. First we have considered only the optical point and the data in the narrow interval 0.01 < 
q 2 < 0.1 GeV 2 (very forward direction). In the logarithmic scale, the differential cross section 
data follows a well defined straight line, so that the slope and the intercept can be directly 
evaluated from the plot and used as start values for the fit with only one exponential term in 
parametrization (13), that is, n = 1. 

2. Once obtained the final fit result, we have enlarged the interval up to g^ax — 0.2 GeV 2 . Due to 
a sudden change of the slope around 0.13 GeV 2 , we have considered two exponential terms (n 
= 2) and used as start values for the parameters in both terms the final values of the previous 
fit. With this procedure the break of the slope is quite well described, as well as all the data in 
the above interval. 

3. Next we have considered g^ax = 0-5 GeV 2 and this case demands n = 3. As start values for the 
third exponential terms we have tested the final values of each of the two terms in the previous 
step and selected the best result (reduced \ 2 closest to 1). 

4. The same procedure has been then applied for = 1.0, 2.0, 3.0, 4.0 GeV 2 and after that, 

= 14.2 GeV 2 . 

Once obtained the best result at 52.8 GeV, the final values of the parameters have been used as 
feed back for the fits at each nearby energy set, namely from 52.8 GeV to 62.5 GeV, from 52.8 GeV 
to 44.7 GeV and then from 44.7 GeV to 30.7, 23.5 and 19.4 GeV, successively 
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TABLE I: Fit results for each data set and statistical information: degrees of freedom (DOF) and reduced chi 
square (x 2 /DOF). For i=l,2,...,6, the parameters a; are dimensionless and bi in GeV -2 . 



\fs (GeV): 


19.4 


23.5 


30.7 


44.7 


52.8 


62.5 


«i 


0.111 


0.506 


0.333 


0.0791 


0.151 


0.179 




±0.011 


±0.029 


±0.019 


±0.0097 


±0.017 


±0.022 


ai 


5.150 


5.734 


5.740 


5.825 


6.8323 


6.7449 




±0.014 


±0.032 


± 0.023 


±0.017 


±0.018 


±0.019 


03 


-8.4565 


-8.4602 


-8.4491 


-8.444 


-7.5697 


-7.6620 




±0.016 


±0.033 


±0.024 


±0.017 


±0.018 


±0.019 


a -i 


-0.985 


-0.00665 


-0.0089 


-0.0099 


-0.00862 


-0.0064 




±0.011 


±0.00085 


±0.0011 


±0.0012 


±0.00082 


±0.0012 


«5 


-0.00165 


-0.000142 


-0.000271 


-0.00037 


-0.00023 


-0.000203 




±0.00010 


±0.000069 


± 0.000090 


±0.00011 


±0.000075 


±0.000093 


C16 


5.187 


3.219 


3.384 


3.564 


1.598 


1.742 




±0.017 


±0.036 


±0.027 


±0.018 


±0.019 


±0.020 


6l 


14.7 


7.82 


9.56 


31.3 


16.4 


14.1 




±1.7 


±0.45 


± 0.60 


±5.8 


±1.8 


±1.7 


b 2 


3.2466 


3.3610 


3.3095 


4.1790 


3.9811 


3.786 




± 0.0045 


±0.0061 


±0.0040 


±0.0085 


±0.0068 


±0.013 


63 


3.7449 


3.4191 


3.4232 


4.4031 


4.0619 


3.852 




±0.0045 


±0.0033 


±0.0083 


±0.0078 


±0.0069 


±0.013 


In 


2.7230 


1.004 


1.150 


1.242 


1.156 


1.079 




±0.0095 


±0.055 


±0.057 


±0.067 


±0.046 


±0.079 


65 


0.5993 


0.373 


0.433 


0.469 


0.426 


0.418 




±0.0084 


±0.046 


±0.036 


±0.035 


±0.036 


±0.046 


be 


4.277 


3.610 


3.8146 


5.121 


5.1138 


4.756 




±0.010 


±0.015 


±0.013 


±0.025 


±0.057 


±0.078 


DOF 


302 


161 


200 


235 


233 


152 


X 2 /DOF 


3.00 


1.59 


1.42 


2.11 


1.74 


1.23 



As described above, we did not impose any kind of constraint in the fit procedure, in special in 
what concerns the number of parameters or possible dependencies of the parameters with the energy. 
The only goal has been to obtain the best fit on statistical grounds. 



C. Results and Comments 



The final values of the free parameters and statistical information on the fit results for each set 
analyzed are displayed in Table I. With the Minuit code, the error-matrix provides the variances and 
covariances associated with each free parameter. This information, together with the statistical errors 
of CTtot and p, are used in the evaluation of uncertainty regions in the differential cross section and all 
the extracted quantities, through standard error propagation procedures (fJ0| . The results together 
with the corresponding experimental data are shown in Fig. 1 in all the q 2 region and at the diffraction 
peak. The contributions to the differential cross sections from the real and imaginary parts of the 
amplitude, 

-j-tt = 7T [Re A] 2 , ^^[Im^f, (17) 
dq z dq z 

are shown in Fig. 2, where we display only the uncertainty regions and the experimental data. 

Before considering the extracted quantities some comments on the fit results presented in this 
section are in order, as follows. 
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From Table I, although the values of each parameter in each one of the six sets analyzed are 
practically of the same order of magnitude, a dependence on the energy cannot be inferred. That is a 
consequence of our fit procedure, developed without any kind of constraint in the free parameters, as 
well as the relatively small energy interval investigated (~ 19 - 60 GeV). In fact, the same undefined 
energy dependence on the free parameters characterizes the results in our previous analyses, with 
both the constrained Q and unconstrained [|| parametrization. In this respect, the new data on the 





FIG. 2: Experimental data and uncertainty regions in the contributions to the differential cross sections from 
the real (gray) and imaginary (black) parts of the amplitude. 
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differential cross section from the LHC (7 - 14 TeV) reaching the large momentum transfer region 
may bring important information for further investigation on possible energy dependencies. 

Although Figure 1 shows that all the experimental data are quite well described, from Table I we 
note that the values of the x 2 /DOF for about 200 DOF correspond to extremely small confidence 
intervals. That is a consequence of the fact that we have considered only statistical errors and not the 
systematic ones. The inclusion of the systematic errors (for example in quadrature with the statistical 
ones) lead to quite good confidence intervals, but we do not think this procedure has statistical 
meaning, as commented at the end of in Sect. III. A. 

From the plots in Figure 2, the contributions to the differential cross section from Re A and Im A 
are in agreement with the standard or generally expected behavior: dominance of the contribution 
from Im A except at the dip region, which is filled up by the contribution from Re A. In all sets 
Im A present a zero (change of sign) at the dip position and Re A two zeros, one at small values of 
the momentum transfer (in accordance with a theorem by A. Martin for even amplitudes (6l| ) and 
another one around q 2 = 2 GeV 2 . We also note that the large uncertainty regions associated with 
Re A, as compared with those from Im A, are consequence of the derivative term in parametrization 
(13) and the variances/covariances involved. Also, from the plots we see that the uncertainty regions 
from Re A are larger at the lower energies, which is explained by our fit procedure that started at 
52.8 GeV, and also the derivative term referred to above. 

At last it should be noted that some characteristics of the fit results here presented and also of 
those in our previous empirical analysis [6J, |9j are consequences of the data ensemble considered, the 
choice of the analytical parametrization as sum of exponential in q 2 and the use of statistical errors 
only and not systematic ones. In this respect three comments are in order as follows. 

(1) All our analyses have been based on the addition of the data at 27.4 GeV to the five ISR sets 
(23.5 - 62.5 GeV). The empirical justification for this procedure has been discussed in detail by Avila 
and Menon in [9J, Sect. 3.3.2. Specifically, by testing different cutoffs in the momentum transfer, 
it has been verified that for q 2 > 3.5 GeV 2 all the above differential cross section data follows a 
power law in q 2 with power A = —7.85 ±0.04 Q, indicating therefore an energy-independent behavior 
of the differential cross section in this kinematic region. On the other hand, in the QCD context, 
it is expected that this region of large momentum transfer in elastic scattering may be accessed by 
perturbative techniques. In this respect the above result favors the triple-gluon exchange picture by 
Donnachie and Landshoff [57J, l58|. which predicts A = —8, but not the constituent interchange model 
(CIM) by Lepage and Brodsky [62|], since in this case it is predicted A = —10. This approach also 
predicts a scaling of the form 



dq 2 ~ S 



which has been verified at lower energies [621 ] . There is also indication of agreement with the data 
at 19.4 GeV and 27.4 GeV, if a particular power law is assumed for f(q 2 /s) and the systematic 
uncertainties of 15% present in both sets are neglected (56j . This small energy dependence in the 
interface between Fermilab and ISR energies corroborates the analysis by Avila and Menon in the 
sense to include the data at 27.4 GeV only at the ISR sets and not at 19.4 GeV (see [a], Sect. 3.3.2, 
for more details). 

(2) Another aspect in our analyses concerns the assumption of a sum of exponential in q 2 , for the 
whole range of momentum transfer investigated. This choice has two foundations as follows. The 
first one is strictly empirical since for the differential cross sections in a logarithm scale all the details 
of the diffraction peak, the dip-bump structure and the large q 2 region can be well described by 
compositions of straight lines, analogous to an exhaustion process (the empirical character of our 
analysis impose no limits in the number of free parameters or exponential terms). The second one, to 
be treated in the next section, concerns the possibility of analytical translation from the <7 2 -space to the 
impact parameter space (Fourier transform) and, more importantly, the analytical evaluation of the 
uncertainties in some extracted quantities, through error propagation from the fit parameters. On the 
other hand, as referred to above, the justification for the inclusion of the data at 27.4 GeV to the ISR 
sets was based on a power law fit analysis at the large momentum transfer region and not exponential 
forms. In this respect, it should be noted that as the momentum transfer increases, a replacement of 
an exponential decrease of the differential cross section to a power one may be related to a transition 
from soft to hard sectors and in the phenomenological context can be connected with the geometrical 
structure of the Regge trajectories, from linear to non-linear one, respectively, as discussed by Fiore 
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et al. [3\ (Sect. 6). Therefore, it might also be important to investigate the possibility to combine 
exponential and power forms in an unique analytical parametrization. However, we shall limit our 
analyses to exponential forms for the reasons referred to above, namely they represent a suitable 
empirical ansatz, "unbiased by any theoretical prejudice" , as suggested by Fiore et al. @ and allows 
analytical translation to the impact parameter space. 

(3) The last aspect concerns the dataset considered and the use of statistical errors only and not 
systematic ones. In this respect it should be noted that another data sets exist, as that compiled and 
analyzed by Cudell. Lengyel and Martynov [63| (where the systematic errors are taken into account) 
and available in [64|. We think it would also be interesting to develop new data reductions with this 
data set and our proposed parametrization; a comparative analysis on all the results could be useful. 



IV. EXTRACTED OVERLAP AND EIKONAL FUNCTIONS 



In this section we consider some implications of the fit results in what concerns properties of 
the extracted overlap functions and the imaginary part of the eikonal (opacity) in the momentum 
transfer space. We first recall the main formulas associated with the impact parameter and eikonal 
representations (Sect. IV. A) and then discuss in some detail the extracted inelastic overlap function 
(Sect. IV. B) and the opacity function in the momentum transfer space (Sect. IV. C). 

A. Impact Parameter and Eikonal Representations 

This subject is treated in detail in In what follows we only recall the main formulas of interest 
here. 

The representation of the elastic scattering amplitude in the impact parameter space is named 
profile function and in case of azimuthal symmetry they are connected by 

/>oo 

A(s,q)=i bdbJ (qb)r{s,b), (18) 
Jo 

where b is the impact parameter, T(s, b) the profile function and Jo is the zero-order Bessel function. 
In the eikonal representation the profile is given by 

T(s,b) = l-e 1 *^, (19) 

where x(s, b) is the eikonal function. 

The unitarity principle in the impact parameter space is usually expressed in terms of the total, 
elastic and inelastic overlap functions, 

G tot (s,b)=G el (s,b) + G in (s,b), (20) 
which, in terms of the profile function reads 



2Rer(s,b) = |r(s,b)| 2 + G in (s,b), (21) 

where G- ln (s, b) is the inelastic overlap function and it follows that the integrated inelastic cross section 
is given by 

/•OO 

a in ( S )=27r/ bdbG m (b,s). (22) 
Jo 

In the eikonal representation, from Eqs. (19) and (21), 

G in (s,6) = l-e- 2Imx ( s ^. (23) 

Since unitarity implies Im x(s,6) > we have G[ n (s,b) <1, so that from (22) G ln can be interpreted 
as the probability of an inelastic event to take place at given b and s; in the black disk limit G- ln — > 1. 
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From the above result, the imaginary part of the eikonal is associated with the absorption in the 
scattering process and for that reason is named opacity function, which we shall denote 

Q(s,b)=Im X (s,b). (24) 

In the momentum transfer space it is given by the symmetrical Fourier-Bcsscl transform, with notation 

p oo 

Cl(s,q) = bdbJ (qb)n(s,b). (25) 
Jo 

Eikonal models can be classified or distinguished according to different choices for the opacity function 
in the momentum transfer space, which is expected to be connected with form factors and elementary 
amplitudes. 

Due to the importance of the physical interpretations associated with the inelastic overlap function 
and the opacity function in the momentum transfer space, we shall display and discuss these extracted 
quantities from the fit results. However, before that, we stress again that the relatively small interval 
investigated (~ 19 - 60 GeV) and the unconstrained fit procedure, turn out difficult to infer energy 
dependencies on these quantities since, in general, the error propagation leads to an overlapping of 
the uncertainty regions. We understand that this does not constitute a serious drawback in our 
analyses because in the phenomenological context what is theoretically unknown, in general, concerns 
the impact parameter and momentum transfer distributions and not the dependencies on the energy 
0-[3[ . Therefore we shall limit the results to typical cases of interest at fixed energy, or those for which 
there is no overlapping of the uncertainties. 



B. Inelastic Overlap Function 



With parametrization (13), the profile function can be analytically evaluated through the inverse 
(symmetrical) transform of Eq. (18), 

I* oo 

T(s,b) = -i qdqJ (qb)A( S ,q), (26) 
Jo 

as well as the overlap functions in Eqs. (20) and (21), together with the associated uncertainties 
through error propagation. For each 6-interval considered for the extracted quantities, we have gen- 
erated 10 3 empirical points in the form of a central value and the corresponding uncertainty, namely 
G ± AG. 

The results for the extracted overlap functions in all the energies investigated are very similar and 
are illustrated in Fig. 3 in the case of y/s = 52.8 GeV (largest set of experimental data). The figure, 
in linear scale, shows the uncertainty regions for the extracted total, elastic and inelastic overlap 
functions. We note that at b = 0, G;„ reaches « 92% of the black disc limit, as expected [65l. l66| . 

We now focus on the extracted inelastic overlap function. In order to develop some test on the 
reliability of the extracted results, we have evaluated the integrated inelastic cross section through 
Eq. (22), a quantity that was not used as experimental input in the data reductions. The results 
for all the sets analyzed are displayed in Table II, showing plenty agreement with the corresponding 
experimental data. The data at 19.4 GeV has been evaluated from er tot = 38.98 ± 0.04 mb obtained 
by Carroll et al. and a e i = 6.78 ± 0.13 mb by Schiz et al. Those at the ISR region come 

from the analysis by Amaldi and Schubert [54l |. 

Next, let us discuss the structure of G; n in the impact parameter space. The result at i/s = 52.8 
GeV, in the form of points with errors and in a logarithm scale, is displayed in Figure 4. First we 
note that in the region - 2 fm the empirical points follow a Gaussian dependence with center at 
the origin, but above this region the decrease is slower and deviate from the Gaussian in the form 
of a tail. This effect is not new since it has already been observed by Amaldi and Schubert [54| and 
also Henyey, Tuan and Kane 69] in the seventies. In terms of the differential cross section data it is 
connected with the change of the slope at « 0.13 GeV 2 , referred to in Sect. III.B. 

In order to go further, looking for some quantitative information on this tail effect, we have tested 
different analytical parametrization for the empirical points in Fig. 4 and in this case we have also 
considered all the energies investigated. Making use of the CERN-Minuit code and using as test 
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TABLE II: Integrated inelastic cross section evaluated through Eq. (22) with the extracted inelastic overlap 
functions and the experimental data [54l. Ir37l. lo8| . 





fit 
< 


exp 


(GeV) 


(mb) 


(mb) 


19.4 


32.03±0.04 


32.11±0.14 


23.5 


32.21±0.16 


32.21±0.14 


30.7 


32.99±0.13 


32.98±0.14 


44.7 


34.62±0.11 


34.62±0.14 


52.8 


35.23±0.13 


35.22±0.16 


62.5 


35.64±0.16 


35.66±0.21 



functions exponential, Gaussian and combinations of these functions, the best global result has been 
obtained with a composition of three Gaussian, which we shall denote 

G in (b) = Yl Ae- B ^ b - C ^ 2 = Gi(b), (27) 

i=l i=l 

where Ai, Bi and Ci, i = 1, 2, 3 are free fit parameters. The result of the fit, through the CERN- 
Minuit code, for pp at 52.8 GeV is shown in Fig. 4, together with the corresponding components 
Gi(b), i = 1, 2, 3; the values of the parameters are displayed in Table III. In this case, since the 
ensemble corresponds to extracted empirical points with errors and not to experimental data we did 
not evaluate the uncertainties associated with these components and make no reference to the goodness 
of fit (reduced % 2 ). 

The results for each Gaussian component, Gj, i = 1, 2, 3, in each energy investigated are displayed 
in Fig. 5. Although, once more, a dependence on the energy cannot be extracted in this interval, we 
can identify the following remarkable features: (1) each component is concentrated in specific regions 
of the impact parameter space, with centers at b = (Gi), around b — 0.6 fm (G2) and roughly 




FIG. 3: Uncertainty regions of the extracted overlap functions for pp scattering at 52.8 GeV and the black 
disc limit (dashed line). 
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FIG. 4: Extracted points with uncertainties for the inelastic overlap function at 52.8 GeV. Left: fit with one 
Gaussian centered at the origin (dashed) and parametrization (27) (solid). Right: result of the fit and the 
three Gaussian components in Eq. (27). 



TABLE III: Fit results for the Gaussian components of the inelastic overlap function, Eq. (27) at 52.8 GeV. 



i 


Ai 


B, 


a 






(fin" 2 ) 


(fin) 


1 


0.9005 ± 0.0067 


0.908 ± 0.020 


0.0 ± 0.0 


2 


0.0866 ± 0.0048 


3.71 ± 0.32 


0.655 ± 0.015 


3 


0.00821 ± 0.00020 


0.684 ± 0.024 


1.653 ± 0.034 




b(fm) MM Mfm) 



FIG. 5: Gaussian components of the inelastic overlap function, Eq. (27), for each set analyzed. Solid lines 
correspond to the extrem energies (19.4 and 62.5 GeV) and the other symbols to the intermediate energies. 
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TABLE IV: Contribution from each Gaussian component to the integrated inelastic cross section and the sum 
of the contributions. 











V^3 Gi 


(GeV) 


(mb) 


(mb) 


(mb) 


(mb) 


19.4 


29.67 


2.87 


0.59 


33.13 


23.5 


29.02 


2.48 


1.29 


32.79 


30.7 


29.23 


3.01 


1.63 


33.87 


44.7 


31.02 


3.52 


1.60 


36.14 


52.8 


31.15 


3.31 


1.84 


36.30 


62.5 


30.26 


4.26 


2.12 


36.49 



around b = 2.0 fm (G3); (2) Gi dominates the central and intermediate regions and G3 describes the 
tail effect above b ss 2.2 fm; (3) the effect of different energies is more evident in G3 then in the others 
components. 

In order to get some quantitative information on the relative contribution of each Gaussian compo- 
nent in the integrated cross section, we have evaluated the following partial quantities: 



<T? i (s) = 2TT bdbG l (b,s), i = 1,2,3. (28) 
Jo 

The results with the corresponding sum of the partial contributions in each energy are shown in Table 
IV. Certainly this sum of the areas under the curves cannot be identified with o~- ln , but allows to 
infer that, roughly, the relative contributions read 85 % from Gi, 10 % from Gi and 5 % from G3. 
Moreover, from Table IV, in this energy interval the relative errors associated with the average of 
each contribution, namely A < aQ > / < >, correspond to 3%, 19% and 35% for i = 1,2,3 
respectively. These quantitative results can be interpreted as evidence that the peripheral region is 
more sensitive to variations of the energy than the intermediate region and this one yet more sensitive 
than the central region. That seems a reasonable result, since we are treating the inelasticity associated 
with the elastic channel through unitarity. Also, from Table IV we can see that except for the result 
at 44.7 GeV, the contribution from G3 increases with the energy. 

In conclusion, these results indicate three Gaussian components for the inelastic overlap function 
associated with central, intermediate and peripheral regions. Although an explicitly dependence on 
the energy cannot be obtained the results suggest that the peripheral region is more sensitive to the 
energy evolution then the central region. We understand that the correct evaluation of <7; n (Table II) 
and the reproduction of the tail effect in Gi n give support to our representation and fit procedure. 



C. Opacity Function in the Momentum Transfer Space 



With the extracted profile function we determine the opacity function in the impact parameter 
space analytically, through Eqs. (19) and (24), 



Q(s, b) = In • 



1 



v / [l-Rer(s,6)] 2 + [Imr(s,fe)] 2 



(29) 



together with the corresponding uncertainties by error propagation from the fit parameters. 

Now, since T(s,b), from A(s,q) in (13), is given by a sum of Gaussian in b, the translation to q 2 - 
space, through Eq. (25), cannot be analytically performed and therefore the error propagation neither. 
For that reason we have used the semi-analytical approach discussed by Avila and Menon Q which 
is valid under certain condition (as explained in what follows) and that we shall name conditioned 
expansion method. 

From the fit results and in all the energies analyzed, we have found that 



Imr(s, b) 
1 -RcT(s,6) 



= r(s,b) < 1, 



(30) 
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including the uncertainty regions. This is illustrated in Fig. 6 in the typical case of y/s = 52.8 GeV. 
Therefore, Eq. (29) for the opacity can be approximated by 



Q(s,b) w In 



1 



1 - ReT 



(31) 



In order to avoid loss of information, we look for what can be directly obtained from the fit. In this 
case, by adding and subtracting Re T, the above equation can be expressed by 



n(s,b) = ReT(s,b)+R(s,b) 
where the "residual" function is given by 

1 



R(s,b) = 



In 



ReT(s,6) 



1-ReT 

Now, the Fourier-Bessel transform of (32) reads 

fi(s, q) = ImA(s, q) + R(s, q) 



(32) 



(33) 



(34) 



so that the first term is directly obtained from the fit and the residual function can be determined 
through a semi-analytical method [9]. Specifically, with (33) we generate a set of points with errors 
(from propagation), Ri(s, b) ± Ai?^(s, 6), i = 1, 2, ..N and parametrize this set by Gaussian forms, 



R{s,b)=Y J A j t 

3=1 



-Bib 2 



(35) 



which has analytical Fourier-Bessel transform and allows error propagation. Then, substitution in 
(34), together with Iim4 from the fit, leads to Cl(s, q) ± Af2(s, q). 

Following the procedure by Avila and Menon [9j we display the results in terms of q s multiplied by 
i7(<7) together with the uncertainty regions, as shown in Fig. 7 for all sets investigated. A remarkblc 



0.0014 



0.0012 



0.001 



0.0008 - 



0.0006 



0.0004 



0.0002 




b(fm) 



FIG. 6: Empirical points at 52.8 GeV for the ratio r(s, b) in Eq. (30) and uncertainties from error propagation. 
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feature in these results is the evidence of a zero (change of sign) in the opacity at finite values of 
the momentum transfer (around 6.5 GeV 2 ), since even the uncertainties in q 8 £l(s,q) lie below the 
zero above the crossing point. Beyond this change of sign, the opacity goes smoothly to zero through 
negative values. 

From the plots in Fig. 7 the position of the zero can be determined and the error involved estimated 
from the uncertainties in both sides of the zero position. The results for these zero positions in terms 
of the energy are displayed in Fig. 8. We note that in the previous analysis by Avila and Menon 
with the constrained parametrization [9| there has been no evidence for a zero at 19.4 GeV. That is 
a consequence of the new compilation and normalization for this data set used here and described in 
detail in |6J. 

In order to obtain an analytical parametrization for the empirical points displayed in Fig. 7, we 

consider a general structure inspired in eikonal models with zero in the opacity function (70l l7l| and 
based on the Glauber multiple diffraction formalism (72l - [74l |. For fixed s we express: 






FIG. 7: Extracted opacity function in the momentum transfer space multiplied by q and uncertainty regions 
from error propagation. 



D.A. Fagundes and M.J. Menon 



19 



with 



and 



ft(g) = CG 2 (g)/(g), (36) 



G(<?) (l + g 2 /a 2 ) (l+g 2 //3 2 ) (3?) 



where C, a, /?, go are free fit parameters and n a fixed integer. The change of sign (zero) in £l(q) is 
generated at q 2 — q 2 . 

These are general and useful formulas, since for different values of the integer n in /(g) they 
reproduce different model assumptions or empirical parametrization for the corresponding opacity 
function. In fact, for n = 1 we have a structure used in the Bourrely, Soffer and Wu model [75T - [78| . 
which we shall denote flssw- The case n = 2 corresponds to a modification in the BSW function, 
which appears in multiple diffraction [H, [49|, [Ho| and hybrid [79| models and we shall denote £l m BSW 
(standing for modified-BSW) . At last, in the analysis by Avila and Menon the value n = 4 has been 
inferred on empirical basis and for that reason we shall denote £t e mpir- 

As before, in what follows we shall focus on the typical results obtained at y/s = 52.8 GeV. The 
point is to fit the empirical points with errors for the opacity function, as displayed in Fig. 7, through 
parametrization (36-38), in the cases n = 1 (Qbsw), n — 2 (Q m BSw) and n = 4 (£l em pir)- In this 
procedure we have fixed the position of the zero to the empirical value at this energy, namely gg = 
6.34 GeV 2 and let free the parameters C, a 2 and (3 2 . The results of the fits are displayed in Table V 
and Fig. 9, for q s (l and also the modulus |f2|. 

We conclude with this analysis that only fl em pir can reproduce the empirical points and in this 
case, from Table 5, /3 2 s» 1 + a 2 (a form factor discussed by G. Shaw, in the beginning of the seventies 
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FIG. 8: Position of the zero in the opacity function extracted from the plots in Fig. 7 (see text) in terms of 
the energy. The dashed lines are drawn only to guide the eyes. 
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TABLE V: Results of the fits to the empirical points in Fig. 7 through Eqs. (36-38) at -0s = 52.8 GeV. 



parameters 


O.BSW (n = 1) 


ClmBSw (n = 2) 




C (GeV" 2 ) 


14.101 ± 0.024 


12.303 ± 0.022 


11.162 ± 0.031 


a 2 (GeV 2 ) 


0.5515 ± 0.0024 


0.6509 ± 0.0040 


0.4622 ± 0.0041 


P 2 (GeV 2 ) 


0.5515 ± 0.0024 


0.6509 ± 0.0040 


1.431 ± 0.018 


X 2 /DOF 


136 


40 


0.19 



[i3,[8l|). These results are in plenty agreement with those obtained by Avila and Menon through the 
constrained parametrization Q and favor models with a zero in the eikonal (see Q for more details) . 



V. CONCLUSIONS AND FINAL REMARKS 



The absence of a pure QCD treatment of the elastic hadron scattering and the lack of a widely 
accepted phenomenological approach to these processes, have motivated us to investigate the inverse 
problem, in the context of the impact parameter and eikonal representations. As commented along 
the paper two of the main problems with these model-independent analyses are: (1) we do not have 
information on the contributions from the real and imaginary parts of the amplitude beyond the 
forward direction; (2) the fits being non-linear have no unique solution. As a consequence, different 
kinds of solutions and contributions from the real and imaginary parts of the amplitude must be tested 
and reliable information should not strongly depend on the details of these solutions and contributions. 




FIG. 9: Extracted opacity function at \/s = 52.8 GeV and parametrization through Eqs. (36-38) for different 
values of n in Eq. (38). Numerical and statistical information are shown in Table V. 
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With this strategy in mind, we have already developed two different empirical solutions for data 
reductions both statistically consistent and characterized by different contributions from the real and 
imaginary parts [H, Q . Motivated by these differences we have discussed here the applicability of a 
representation for the Martin's formula without the full scaling property, which predicts also different 
contributions from the previous analyses. 

A detailed review on the original derivation of this formula and critical comments on its further 
developments have been presented. The use of the proposed almost model-independent representation 
in empirical fits of the differential cross section data in the interval 19.4 GeV - 62.5 GeV has led to 
the following main conclusions: 

• In the interval investigated, the experimental data on the differential cross sections at q 2 = 
and q 2 > 0.01 GeV 2 are quite well described. 

• The extracted inelastic overlap function cannot be parametrized by one Gaussian distribution 
due to the tail effect, associated with the change of the slope in the differential cross section at 
small values of the momentum transfer. The empirical results indicate three Gaussian compo- 
nents whose sensitivity to variations of the energy increases from the central to the peripheral 
regions. 

• The opacity function in the momentum transfer space, extracted through the conditioned expan- 
sion method, present a change of sign around 6.5 GeV 2 in all sets analyzed. The empirical points 
are quite well described by a parametrization consisting of a product of two simple constrained 
poles (/3 2 w 1 + a 2 ) by the function f(q) with zero and n=4. 

Despite the promising aspects of the above results we stress that only a detailed comparative 
investigation, based on different empirical parametrization and solutions and aimed to identify results 
that are common to all the analyses, can bring new reliable insights for model developments and 
possible connections with QCD. A global critical review along these lines, collecting all our previous 
results and those presented here, are being prepared for a forthcoming communication. 
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